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Abstract. We employ a long XMM-Newton observation of the core of the Perseus cluster to validate claims of a non- 
thermal component discovered with Chandra. From a meticulous analysis of our dataset, which includes a detailed treat- 
ment of systematic errors, we find the 2-10 keV surface brightness of the non-thermal component to be smaller than about 
5x 10^'^erg cm^-s"' arcsec"^. The most likely explanation for the discrepancy between the XMM-Newton and Chandra estimates 
is a problem in the effective area calibration of the latter. Our EPIC based magnetic field lower limits are not in disagreement 
with Faraday rotation measure estimates on a few cool cores and with a minimum energy estimate on Perseus. In the not too 
distant future Simhol-X may allow detection of non-thermal components with intensities more than 10 times smaller than those 
that can be measured with EPIC; nonetheless even the exquisite sensitivity within reach for Simbol-X might be insufficient to 
detect the IC emission from Perseus. 
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1. Introduction 

Although the bulk of the energy radiated from clusters is of 
thermal nature, non-thermal mechanisms play an important 
role. Indeed the characterization of non-thermal components 
provides much needed clues on the physical process presid- 
ing over the formation and evolution of clusters. In some 
of the more disturbed objects, evidence of non-thermal pro- 
cesses has been known for quite some time. Radio observa- 
tions indicate that merging clusters are often the site of cluster- 
wide sync hrotron emission, th e so called radio halos and ra- 
dio relics dCassano et al.ll2007l and references therein). Radio 
data, polarimetric or not, has been used to provide estimates of 
magnetic fields in clusters. Unfortunately estimates based on 
Faraday rotation measures (hereafter RM) or minimum-energy 
arguments are affected by large uncertainti es. In the case of RM 
estimates the unknown field topology (.EnBUn & Vogj |2003|) 
and the accessability of only a few, not necessa rily represen- 
tative, lines of sight ( Rudnick & Blundellll2003h . are a major 
source of concern, while for minimum energy arguments the 
proton to electron ratio, and the applicability of the argument 
itself play an equally critical role. 

It has been recognized for quite some time that detection of 
inverse Compton (hereafter IC) emission at X-ray wavelengths 
can pr ovide an alternative m ethod to estimate cluster magnetic 
fields (IRephaeli et al.lll987b . Detection and characterization of 
the so called hard tails is by no means a trivial task: the signal. 



if there is one, is caught between the hammer of the thermal 

emission and the anvil of the instrumental background. So far, 

only detections at a few sigm a level h ave been reported (for a 

recent review see Rephaeli et al.l2008l) . and at least in one case, 
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Coma (Fusco-Femiano et al. 1999, 2004, 2007), they have been 

challenged (.Rossetti & MolendL2004..2007} . 

Aim of the present work is to employ a long XMM-Newton 
observation of the Perseus core to validate claims about a 
non-thermal emission component discovered wit h Chandra 
(Sanders et al. 2004, 2005; Sander s & Fabianl2007h . The major 
advantage of the XMM-Newton observation is the better sensi- 
tivity at high energies (~ 6 keV) of the EPIC detectors with 
respect to the ACIS-S3\ this is illustrated in Fig. [1] where we 
show the ratio of background to source intensity in a represen- 
tative region, a ring with bounding radii l'-2' centered on the 
emission peak. The Chandra ratio rises very rapidly from 1% 
at 5keV to unity at 8 keV; at 10 keV the ACIS-S3 background 
is about 20 times larger than the source. MOS and pn ratios 
rise above 1% around 7 keV and remain below 30% and 10% 
respectively at all energies. 

The outline of the paper is as follows: in Sect.|2]we provide 
details on our data preparation; in Sect. [3] we describe results 
from the spectra analysis; in Sect. |4] we compare results from 
the analysis of EPIC data with those from Chandra; in Sect.|5] 
we discuss our results and in Sect. |6] we summarize them. We 
assume Hq = 70 km s^'Mpc"' so that 1 arcsec=0.35 kpc at the 
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Fig. 1. Ratio of background to source spectra for the ACIS-S3 
(blue), M0S2 (green) and pn (red) instruments from the r-2' 
ring. The ACIS-S3 ratio rises rapidly from 1 % at 5keV to unity 
at 8 keV; at 10 keV the ACIS-S3 background is about 20 times 
larger than the source. MOS and pn ratios rise above 1 % around 
7 keV and remain below 30% and 10% respectively at all ener- 
gies. 



redshift of NGC1275 (0.0176. iHuchra et al.lll999h : all errors 
are Icr unless otherwise stated. 

2. Data preparation 

Perseus was observed by XMM-Newton during revolution 1 125 
with the EPIC MOS {pn) detectors in Full Frame Mode 
(Extended Full Frame) with medium filter on all EPIC cameras. 
We have obtained calibrated event files for the MOSl, M0S2 
and pn cameras using SAS version 7.1 'Emchain", "Epchain" 
tasks and calibration files (CCFs) available at the end of August 
2007. A soft proton cleaning procedure was performed by ex- 
tracting a light curve in 100 second bins in the 10-12 keV 
energy band in a peripheral region of the MOS and pn cam- 
eras and applying a threshold of 0.20 cts s ' for MOS and 
0.6 cts s"' for pn, to generate a filtered event list. The total 
live time after screening is 1.12 x 10^ s for MOSl and M0S2 
and 6.6 x 10"* for pn. Finally, we filtered event files accord- 
ing to FLAG (FLAG==0) and PATTERN (PATTERN<4 forpn and 
PATTERN==0 for MOS). The filtered event files, which were 
subsequently used for data analysis, contained a total of about 
16 X 10^ 19 X 10* and 35 X 10^ events for MOSl, M0S2 and 
pn respectively. 

2.1. Spectra preparation 

We have accumulated spectra from annuli with bounding radii 
0.5'-l', r-2' and 2'-3', and from sectors 45 degrees wide from 
the same annuli. 

As already pointed out in Sect. |2] for the MOS detec- 
tors only PATTERN==0 events were selected, this was done 
to avoid pile up which, albeit small (few %), can have non- 
negligible consequences within ~ 2' from the cluster core (i.e. 
be sufficiently str ong to distort the spectral shape, particularly 
at high energies, Molendi & Sembavll2003[) . We generate flux 
weighted eff'ective areas using the arf gen task with exposure 



corrected images of the source and the parameter extended 
source switched to true. We also correct pn spectra for ou t 
of time events following the prescriptions of iGrupd (1200 lb . 
Redistribution matrices were generated with the rmfgen task, 
one matrix was generated for MOSl and M0S2 respectively, a 
set of 10 matrices were generated fox pn to account for the de- 
pendence of the pn spectral resolution upon the distance to the 



Although Perseus is the brightest cluster in the X-ray sky 
and EPIC offers the largest effective area of any X-ray exper- 
iment flown thus far, a reliable characterization of the back- 
ground is one of the fundamental requirements to carry out a 
search for non-thermal emission. Our background subtraction 
procedure makes use of a sum o f blank fields selected by our 
own group (ILeccardi & Molendill2008l) . Spectra have been ac- 
cumulated from the same sky regions as the source spectra, af- 
ter re-projection onto the sky attitude of the source. To allow 
for differences in the intensity of the background components 
between source and background observations we re-normalize 
the background spectra. The re-normalization constant is de- 
termined by comparing the intensity of the total emission in 
source and background observations from a peripheral region 
and in a hard (10-12 keV for MOSl and M0S2, 12-14 keV for 
pn) energy band where source emission is expected to be small. 
The re-normalization constants, Rk, are found to be 1.64, 1.57 
and 1.52 for MOSl, M0S2 and pn respectively. To verify the 
correctness of our procedure we performed a spectral fit to the 
data in the outer ring used to determine the re-normalization 
constants, the fit was performed in a hard, albeit source dom- 
inated band (4-7 keV). The best fitting model was then ex- 
trapolated to a background dominated band (8-14 for pn and 
8-12 for MOS) where it was found to reproduce adequately 
the source spectrum. Variations of the re-normalization con- 
stant described above of only a few percent result in a best- 
fitting model extrapolation significantly under or over-shooting 
the source spectrum, indicating that our determination of the 
re-normalization constants is good to a few percent. To ac- 
count for possible variations in the background intensity over 
the detectors and the issuing re-normalization errors on the 
background subtraction we have assumed a 5% systematic un- 
certainty on the re-normalization constants we apply on all 
our source spectra. In practice this is achieved by perform- 
ing 3 different fits for each spectrum one with the background 
with the standard re-normalization constant, Rk, the other two 
with a larger, Rk + 0.05 x Rk, and smaUer, Rk - 0.05 x Rk, re- 
normalization constants respectively. 



As a further check of the quality of our background sub- 
traction technique we have also pursued a different venue , i.e. 
background modeling, following the l ines described in a re- 
cent paper (ILeccardi & Molendil 120081) . A comparison of re- 
sults obtained with the two methods shows that within some 4 
arc-minutes of the cluster core they provide virtually indistin- 
guishable results. Given that the background modeling alterna- 
tive is more cumbersome from the point of view of the spectral 
analysi s and that it has only been developed for the MOS de- 
tectors ( ILeccardi & Molendill2008h in the following we report 
results obtained with the background subtraction method. 
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One of the hallmarks of the Perseus observations is its high 
statical quality, typical spectra contain in the order of hundred 
of thousands to millions of source counts. Under these ex- 
treme conditions statistical errors can be minuscule, especially 
around 1 keV where effective areas peak: here systematic errors 
become an important if not dominant component of the error 
budget. This can be readily verified by inspecting spectral fits, 
typically residuals appear at the position of edges where rapid 
variations in the effective areas occur Conversely far less statis- 
tics is available in the 5-10 keV band where some astrophysi- 
cal features, such as the thermal exponential cutoff, Fe and Ni 
emission lines are located. The question is then how to pre- 
vent the x^ fitting procedure from giving more importance to 
statistically significant spurious features located around 1 keV 
than to astrophysical signatures at higher energies. We adopt 
two independent methods: the first is to identify the cause of 
systematic errors and evaluate their impact on parameter es- 
timates; the second is to associate a systematic error to spec- 
tra. While the former strategy may be considered as part of 
the modeling process of Perseus spectra, and as such will be 
described in detail in the following section, the latter is accom- 
plished by associating a systematic error to all spectral bins. We 
choose a somewhat moderate and energy independent value of 
2% for the systematic error for three reasons: firstly it is consis- 
tent with a conservative estimate for the relative calibration of 



the EPIC instruments dKirsch & EPIC-Consrtiumll2007l) : sec- 
ondly it has the effect of increasing errors in the low energy 
part of the spectrum to the point that they become comparable 
to errors at higher energies; thirdly a detailed energy depen- 
dent characterization of systematic errors is in part achieved 
through the modeling processes described in the following sec- 
tion and in part by associating the same relative systematic er- 
ror to spectral regions characterized by very different statistical 
errors as already clarified in the previous point. Inclusion of 
the systematic error is technically achieved by modifying the 
SYS_ERR keyword in the source spectra header. The XSPEC 
package associates to each spectral bin an error which is the 
sum in quadrature of the statistical and systematic error. 

3. Spectral analysis 

Due to its higher effective area we have considered the 
pn camera as lead instrument in our analysis. Results from 
the analysis of MOS data is also reported (see Sect. 13.1.2b . 
We have made use of two kind of baseline models: a 
two temperature plus power-law model, in XSPEC jar- 
gon phabs'''(2vmekal+pegpwrlw), which we shall refer to 
as 2TP, and a multi temperature plus power-law model, 
phabs* (Svmekal+pegpwrlw), hereafter MTP For the former 
model the free parameters are the Nh, the temperatures, the 
normalizations of the thermal and pow-law components, the 
redshift, the O, Mg, Si, S, Ar, Ca, Fe and Ni abundances. Ne, 
Na and Al abundances are tied to O; metal abundances for a 
given species of one component are tied to those of the other; 
the photon index of the power-law is fixed to 1.8, one of the 
valu es adopted in ISanders et al.l (120051) and ISanders & FabianI 
( 12007.) . For the latter model the free parameters are the Nh, the 
normalizations of the thermal and power-law components, the 




Energy [keV] 

Fig. 2. Best fit with the 2TP model to the pn spectrum in the 
l'-2' ring. In the top panel we show the spectrum, in black, 
and the best fitting model in red. In the bottom panel we show 
residuals in the form of a ratio of the data over the model. 



redshift, the O, Mg, Si, S, Ar, Ca, Fe and Ni abundances. Ne, 
Na and Al abundances are tied to O; metal abundances for a 
given species of one component are tied to those of the oth- 
ers. The five temperatures, which are fixed, are kTi = 0.5 keV, 
kT2 = 1.0 keV, kTa = 2.0 keV, kT4 = 4.0 keV and kTj = 8.0 
keV, the power-law index, as for the 2TP model, is fixed to 
1.8. Throughout our analy sis we adopt solar units according to 
Anders & Grevessd (Il989h . 

Spectral fits are conducted in the 0.5-1 1 keV band for both 
MOS and pn. Th e low energy cutoff is enfor ced to avoid resid- 
ual calibrat ion ( Kirsch & EPIC-ConsrtiumI l2007h and cross- 
calibration (IStuhlinger et al.l2008b problems, moreover the rel- 
atively high column density towards the Perseus cluster reduces 
significantly the value of data accumulated at the softest ener- 
gies. 



3.1. Analysis in rings 

We investigate a region within 3' from the emission peak, 
where evi dence of non-thermal em ission has been found with 
Chandra (ISanders & Fabianll2007l and references therein). We 
consider spectra in 3 annuli with bounding radii: 0.5'- 1.0', 1.0'- 
2.0' and 2.0'-3.0'. We exclude the region within 0.5' from the 
emission peak to avoid contamination from the AGN. 



From the analysis of the Chandra data ( ISanders & Fabian 



20071 and references therein) we expect the putative non- 
thermal component to make up no more than 10-20% of the 
total emission, even at high energies where its relative contri- 
bution is largest. Under such conditions a meticulous analysis 
of the data is mandatory. We begin by considering the annulus 
with bounding radii r-2'. 

3.1.1. The 1 '-2' ring 

A fit to the pn data with the MTP model returns values of 
the surface brightness of the non-thermal component, Snt, of 
(1.9+0.3)x 10^'^erg cm^^s^'arcsec"^, the 2TP model returns a 
similar value of (1.6 ±0.3) xlO~'Srgcm"^s"'arcsec"^ (see also 
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Table 1. Parameters from the pn l'-2' spectrum fits 



Model 



2TP-nosys 
MTP-nosys 



2TP 
MTP 



2TP-gs 
MTP-gs 



2TP-gs-2Fe 
MTP-gs-2Fe 



Fe 



0.53: 
0.53: 



:0.01 
:0.01 



0.48: 
0.50: 



0.01 
0.01 



0.50: 
0.52: 



0.01 
0.01 



0.63: 
0.79: 



:0.06 

0.10 



7" 



0.47: 
0.49: 



0.02 
0.01 



4.3: 
3.7: 



0.3 
0.5 



1.6: 
1.9: 



0.3 
0.3 



1.5 ±0.8 
1.7 + 0.7 



1.1: 
1.3: 



0.8 

:0.7 



1540 
1489 



636 
631 



453 
449 



447 
437 



dof 



509 
508 



509 
508 



509 
508 



is possible to have two largely independent measurements of 
the Fe abundance, one for the hot phase, contributing mainly 
to the K-shell line emission, and the other for the cool phase 
dominating the L-shell blend emission. We recall that, for most 
cool core clusters, spatially resolved analysis shows that the 
Fe abundance anti-correlates with temperature (i.e. Fig. 1 1 of 
Sanders et alJ l2004 for the Perseus cluster). We have modified 



508 
507 



Not es: ° F e abundance in solar units according to lAnders & Grevessd 
1 119890 . * Surface brightness in units of 10^'^erg cm^^s^'arcsec"^ 



in the 2-10 keV band. 



Tab. [U. Interestingly a fit of the MTP and 2TP models on the 
pn spectrum with no systematic errors applied, hereafter MTP- 
nosys and 2TP-nosys, returns larger values for the non-thermal 
component surface brightness (see Tab.lTJ. 

Inspections of the best-fit residuals reveals three major fea- 
tures, namely: an excess around 1.2 keV; a substantial struc- 
ture around the Fe Ka line and an excess around 8.3 keV (see 
Fig.|2]). The 1.2 keV excess is located in the blue wing of the Fe 
L-shell blend and is most likely due to an inaccurate modeling 
of the L-shell transitions within the mekal code. Indeed fitting 
with the apec rather than the mekal code results in a substan- 
tial reduction of the residuals, with no significant modification 
of the best fit parameters. The excess around 8.3 keV is due 
to background subtraction inaccuracies. We remind our readers 
that in the pn camera the region around 8 keV is populated by 
fluorescence lines of instrumental origin which are not easily 
subtracted. F ortunately the co ntamination in the V -2' region 
is only mild JEhle et al.ll2007l) and the ensuing systematics are 
relatively modest. A further mitigating factor is the relatively 
large statistical error associated to the source spectrum above 
~ 7keV. The structure around the Fe Ka line is due to an incor- 
rect modeling of the pn spectral resolution within the redistri- 
bution m atrix. From the analysis of pn spectra from a sample of 
clusters (iDe Grandi & Molendil l2008. in preparation) we have 
found that the above resolution miss-calibration can be com- 
pensated for by including a multiplicative component that per- 
forms a gaussian smoothing of the spectral model (gsmooth in 
XSPEC). We set the width of the gaussian kernel to be 4 eV 
(FWHM) at 6 keV and assume a power-law dependency of the 
width on the energy with an index of -1. Both the MTP and 
2TP models modified with the gsmooth component, hereafter 
MTP-gs and 2TP-gs respectively, provide a substantially bet- 
ter fit than their non-modified counterparts {t^x^ ~ 200). The 
surface brightness of the non-thermal component is essentially 
unchanged with respect to the previous fits both for the MTP- 



gs, 5 ivr = (1.7 + 0.7) X 10 



""erg cm ^s 



arcsec and 2TP-gs 



model S NT = (1-5 ± 0.8) x lO^'^erg cm 's 'arcsec 

A further modification concerns the treatment of the Fe 
abundance. When fitting multi-phase models it is customary 
to tie the abundances of the different phases; this solution is 
adopted because in most cases it is impossible to derive in- 
dependent metal abundance measures. In the case at hand it 



our two test models to allow for different Fe abundances in 
4he cool and the hot phase. For the 2TP-gs model we simply 
decouple the Fe abundance of the two components while for 
the MTP-gs model we define two Fe abundances: one for the 
cool i.e. 0.5, 1.0 and 2.0 keV components and the other for 
the hot, i.e. 4.0 and 8.0 keV components. As expected for both 
models, hereafter 2TP-gs-2Fe and MTP-gs-2Fe, the Fe abun- 
dance for the cool phase or phases is higher (Zf^ = 0.63 + 0.06 
and Zfc = 0.79 + 0.10 respectively) than for the hot phase 
(Zpe = 0.47 + 0.02 and Z^, = 0.49 + 0.01 respectively). 
Both the MTP-gs-2Fe and 2TP-gs-2Fe models provide a bet- 
ter fit than their single Fe abundance counterparts although 
the improvement is not a dramatic one {Ax^ ~ 10). The sur- 
face brightness of the power-law component is found to be 
Snt = (1.1 + 0.8) X 10^'^erg cuT^s-^ arcsec-^ for the 2TP-gs- 
2Fe model and S nt = (1-3 + 0.7) x lO^^^erg cm^^s"' arcsec"^ 
for the MTP-gs-2Fe model. In retrospect it is fairly easy to un- 
derstand why the decoupling of hot and warm phase Fe abun- 
dances entails somewhat lower intensities for the non-thermal 
component. A high temperature component (kT ~ 8 keV) over 
a small energy range is rather similar to a power-law compo- 
nent, the major discriminant being of-course the line emission. 
If the high temperature component is forced to have a metal 
abundance that is larger than the true one it will not be able to 
reproduce the continuum correctly, such an excess continuum 
will be well modeled by a power-law component. 

We may summarize our results by stating that while the ad- 
dition of systematic errors to the spectrum results in a reduction 
of the intensity of the power-law component of a factor of 2, 
the introduction of some modifications to our base-line models 
2TP and MTP, either to improve the description of the astro- 
physical scenario or to account for residual calibration prob- 
lems, does not result in substantial modifications of the inten- 
sity of the non-thermal component. 

To asses how different choices for the power-law slope 
impact on the estimate of the surface brightness of the non- 
thermal component we have rerun fits with slopes of F = 1.5 
and r = 2.0. We find that surface brightnesses are some- 
what smaller (larger) for F = 1.5 (F - 2.0). For instance, 
in the case of the MTP-gs-2Fe model, the surface bright- 
ness is Snt = (1-2 ± 0.8) x 10^'^erg cm^^s^Wcsec"^ and 
Snt - (2.0 + 0.8) x 10^'^erg cm^^s"' arcsec"^ respectively for 
the r = 1 .5 and F = 2.0 choices. To gauge the effect of system- 
atic errors on our fitting procedure we have experimented with 
values ranging from 0% to 4%. For the MTP model we find that 
the surface brightness of the non-thermal component steadily 
decreases from (3.7 + 0.5) x 10^'^erg cm^^s"' arcsec"^ for no 
systematic error to (1.2 + 0.3) x 10^'^erg cm^^s"' arcsec"^ for a 
4% systematic error To verify whether spatial resolution plays 
a role in the determination of the power-law surface brightness 
we have analyzed spectra from circular sectors. For the 1-2 ar- 
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Table 2. Parameters from the MOSl V-T spectrum fits 



Table 3. Parameters from the M0S2 V-2' spectrum fits 



Model 


7" 


7" 




T 


dof Model 


7" 


^. 


cl, 


x' 


dof 


2TP-nosys 
MTP-nosys 


0.66 + 0.02 
0.63 + 0.01 




5.4 ± 1.0 
4.2 ± 0.5 


951 
1009 


553 2TP-nosys 
552 MTP-nosys 


0.63 ±0.01 
0.56 ±0.01 




7.4 ±0.4 
2.6 ±0.5 


1066 

1228 


555 
554 


2TP 
MTP 


0.62 + 0.02 
0.59 + 0.01 





3.8 ± 1.1 
3.0 ±0.9 


593 
599 


553 2TP 
552 MTP 


0.61 ±0.02 
0.55 ±0.01 


— 


4.2 ±1.5 
3.4 ±0.7 


635 
657 


555 
554 


2TP-2Fe 
MTP-2Fe 


0.67 + 0.04 
0.66 + 0.08 


0.53 + 0.05 
0.58 ± 0.02 


2.7 ± 1.4 
3.0 ±0.9 


591 
598 


552 2TP-2Fe 
551 MTP-2Fe 


0.65 ± 0.04 
0.71 ±0.07 


0.37 ±0.10 
0.52 ± 0.03 


2.8 ±1.5 

2.9 ±0.8 


632 
655 


554 
553 


2TP-2Fe-gf 
MTP-2Fe-gf 


0.67 + 0.05 
0.71+0.12 


0.54 ± 0.05 
0.58 ± 0.02 


2.7 ± 1.5 
2.4 ±1.3 


611 
592 


551 2TP-2Fe-gf 
550 MTP-2Fe-gf 


0.74 ±0.13 
0.80 ±0.15 


0.49 ± 0.07 
0.52 ±0.02 


3.6 ±1.7 
2.9 ±1.5 


611 
615 


553 
552 



Not es: " F e abundance in solar units according to lAnders & Grevessd 
1 119890 . * Surface brightness in units of 10^'^erg cm^^s^'arcsec"^ 



in the 2-10 keV band. 



Not es: ° F e abundance in solar units according to lAnders & Grevessd 
( Il989l) . * Surface brightness in units of 10" 
in the 2-10 keV band. 



'^erg cm ^s 



cmin annulus we consider 8 equally spaced sectors each 45 de- 
grees wide. Analysis of the sectors provides results similar to 
those found in rings albeit with a somewhat reduced statistics. 



unlike for the MOSl and pn detectors, inclusion of a system- 
atic error results in a reduction of the non-thermal component 
surface brightness only for values of the error larger than 3%. 



3.1 .2. The MOS spectra for the 1'-2' ring 

As anticipated we have used MOSl and M0S2 spectra to verify 
pn results. A fit with the simple MTP and 2TP models returns 
values of the surface brightness of the non-thermal component 
between 3 and 4 xlO""'erg cm"^s"'arcsec"^, (see Tabs.|2]and 
[3] for details). These values are in agreement with each other 
and somewhat in excess of those obtained with the pn. Fits 
on the MOSl and M0S2 spectra with no systematic errors ap- 
plied, dubbed 2TP-nosys and MTP-nosys, resulted in a higher 
value of the;^'^; the surface brightness of the non-thermal com- 
ponent is also generally higher, however this not always the 
case (see Tabs. |2] and [3] for details). Decoupling of the hot 
and cool phase Fe abundances results in changes similar to 
those observed for the pn: the x^ suffers a modest decrease; 
the hot phase Fe abundance is found to be smaller than the 
cool phase Fe abundance and the surface brightness of the non- 
thermal component does not vary substantially. Inspection of 
the residuals shows that the MOSl and M0S2 Fe Ka lines are 
not broader than the one included in the best fitting model, con- 
sequently we refrain from applying a smoothing component 
to the spectral model. Conversely we find evidence of some 
residual non-linearity in the energy scale calibration, particu- 
larly in the M0S2 spectrum. We have compensated for this ef- 
fect by allowing for a constant shift in the energy scale, this is 
achieved within XSPEC by performing a gain fit on the best 
fit model, fixing the linear parameter to 1 and leaving the offset 
parameter free. From the resulting best fits, which we refer to 
as 2TP-2Fe-gf and MTP-2Fe-gf, we detect a shift of about 2 eV 
for MOSl and 6 e V for M0S2, other best fit parameters are not 
substantially different from those measured from the previous 
set of fits (see Tabs.|2]and |3]for details). 

As for pn spectra we have experimented with different val- 
ues for the power-law slopes, the systematic eiTor and per- 
formed analysis of spectra in sectors. Our findings are similar 
to those reported for the pn spectra in as far as the power-law 
slope and analysis in sectors is concerned. There are some dif- 
ferences regarding systematic errors. For the M0S2 detector. 



3.1 .3. Comparison between pn and MOS estimates 

Comparing results for MOSl and M0S2 with those from 
pn we find that the surface brightness of the non-thermal 
component gauged with the former experiments is somewhat 
larger than that measured with the latter, however, for most 
measures, the differences are comparable to the errors asso- 
ciated to individual measurements. More specifically, if we 
exclude measurements without correction for systematic er- 
rors we find that pn measures are roughly contained between 
lxlO""'erg cm^^s^'arcsec"^ and 2xl0"'^erg cm^^s^'arcsec"^ 
and MOS measures between 2xl0"'^erg cm^^s^'arcsec"- and 
4xl0"'''erg cm^^s^'arcsec"^. 

If we consider measurements with no correction for sys- 
tematics (i.e. MTP-nosys bst fits) we find that typical errors, 
O.SxlO^'^erg cm^^s^'arcsec"^, are substantially smaller than 
the spread of values described above, indicating that system- 
atic uncertainties afford the dominant contribution to the er- 
ror budget. Under these circumstances providing a statistically 
meaningful confidence range for our measurement of the sur- 
face brightness of the non-thermal component is rather diffi- 
cult, particularly since the approach we have adopted to incor- 
porate systematic errors is of a heuristic nature. We therefore 
resort to an alternative and less ambitious approach: we make 
use of the spread of values measured for the different mod- 
els and different experiments to provide a loose range for the 
surface brightness of the non-thermal component, namely 0- 
SxlO^'^erg cm^^s^'arcsec"^. It is perhaps ironic that a meticu- 
lous analysis of the data, such as the one reported here, has lead 
to an uncommonly loose estimate of the non-thermal compo- 
nent, however, as long as we do not dispose of statistical tools 
allowing us to incorporate systematic uncertainties in analysis 
such as these, our options are rather limited. 

The results we report here are at variance with those that 
we presented at a conference in Garching in August 2006, es- 
sentially 4 factors have lead to these differences. The factors 
are: a new release of the MOS quantum efficiency curves in 
August 2007 leading to a substantial improvement of the MOS 
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Table 4. Parameters from the Chandra ACIS-S3 l'-2' spectrum 
fits 
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Fig. 3. Radial profile for the surface brightness of the non- 
thermal component for the three EPIC cameras: pn is red, 
MOSl is blue and M0S2 is green. Radii have been slightly off"- 
set to improve readability. 



vs. pn cross-calibration and to a revision of MOS spectral anal- 
ysis results; use of PATTERN==0 events only for MOS spectra; a 
more accurate treatment of the backg round foll owing some re- 
cent work within our group (Leccardi & Molendi,2007, ,2008); 
a novel treatment for systematic effects. 



3.2. Radial profile 

In Fig.[3]we show the radial profile for the surface brightness of 
the non-thermal component. The profile is based on the analy- 
sis of spectra accumulated in annuli with bounding radii: 0.5'- 
1.0', 1.0'-1.5', 1.5'-2.0' and 2.0'-3.0', note that the 1.0'-2.0' 
annulus has been split in two to improve the resolution of the 
profile. The reported surface brightnesses have been estimated 
using the MTP-gs-2Fe model for pn and MTP-2Fe-gf model 
for MOSl and M0S2. 



All 



measurements 



5x10 '^erg cm ^s 'arcsec ^ 



lie between and 
For the innermost ring the 
EPIC mean value for the surface brightness of the non-thermal 
component is (2.4 ± 1.0) x 10^'^erg cm^^s"' arcsec"^, 
while for the 1.5'-2.0' annulus we find Snt = 
(2.3 + 0.9) X 10" ""erg cm^^s^'arcsec"^. 

We have verified that similar measurements are obtained 
when fitting spectra with slightly different spectral models, 
i.e. all pn measurements obtained with the MTP-gs model are 
within Icr of those reported in Fig. [3] similarly all MOSl and 
M0S2 measures obtained with the MTP-2Fe model are within 
Icr of those reported in Fig. [31 

Assuming a rather conservative approach, we might say 
that the surface brightness of the non-thermal component is 
somewhere between and 5xlO"'*erg cm"^s~Wcsec~^. This 
is in agreement with what we found from the analysis of the 
1.0'-2.0' annulus (see Sect. I3.1.3I I and may therefore be con- 
sidered as a general result of our analysis. 
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[Notes:] ° Fe a b undan ce in solar units according to 
lAnders & Grevessd ( Il989h . * Surface brightness in units of 
10^'^erg cm^^s"' arcsec"^ in the 2-10 keV band. 



4. Comparison with Chandra 

In an attempt to understand the origin of the difference be- 
tween our measurements and those based on Chandra data, we 
have reduced and a nalyzed the longest AC/S'-^ observation (see 
Fabian et al.ll2006i Tab. 1). The data were reduced with the X- 
ray analysis packages CIAO 3.4 and Heasoft 6.4 in conjunction 
with the Chandra calibration database (Caldb) version 3.4.2 
(using gainfile acisD2000-01-29gainj:tiN0006.fits). To ensure 
the most up-to-date calibration, all data were reprocessed from 
the "level 1" events files, following the standard Chandra data- 
reduction thread^. We applied the standard corrections to ac- 
count for a time-dependent drift in the detector gain and charge 
transfer inefficiency, as implemented in the CIAO tools. From 
low surface brightness regions of the active chips we extracted 
a light-curve (5.0-10.0 keV) to identify and excise periods of 
enhanced background. The spectra were re-binned to ensure a 
S/N of at least 3 and a minimum 20 counts per spectral bin. 
Spectral fits, using both 2TP and MTP models provide val- 
ues consistent with those derived bv lSanders & FabianI (120071) 
(see their Fig. 20) and well in excess of those measured with 
EPIC (see Tab. 4). We have included a 2% systematic error 
to the Chandra spectrum, as done for the EPIC spectra, with- 
out noticeable effects on the non-thermal component, or any 
other spectral parameter. The background file represents well 
the background in the actual observation, as gauged by the 9- 
1 3 ke V band count-rate and even if we experiment changing the 
normalization by +5% the flux in the power-law component is 
comparable with the one measured with the background at the 
nominal value (7.6 ± 0.7 and 8.9+Qg, in the units of Tab. 4 and 
using the MTP model, in the two cases). We investigated possi- 
ble variations in the set of Perseus observations reducing other 
long observations, namely OBS-ID 4948 and OBS-ID 4289, 
in particular the latter provides an important double check be- 
cause it was performed at an earlier epoch and with a different 
level of low-energy degradation. We did not find any difference 
within the errors for the value of the spectral parameters and in 
particular of the non-thermal component between the different 
observations. 

In an attempt to understand the difference between 
Chandra and XMM-Newton measurements, we have com- 
pared the EPIC-pn spectrum accumulated in the r-2' annu- 
lus with the Chandra spectrum extracted from the same region. 



http ://cxc . harvard, edu/ciao/threads/index. html 
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5. Discussion 



From a general perspective our analysis may be viewed as an 
attempt to characterize the spectrum of an X-ray source beyond 
a simple one/two component model. Given the modest intensity 
of the component we are specifically interested in and the lack 
of spectral signatures, such as lines or edges, the task is by no 
means a trivial one. Clearly it can only be attempted on sources 
where adequate statistics is available. It is often the case that, 
when the statistics is abundant, systematics become the dom- 
inant source of eiTor. At the present time X-ray astronomers 
do not have a standard method of analyzing data under these 
predicaments and mostly resort, has we have done, to a trial 
and eiTor approach. The development of a standard strategy, 
possibly descending from first principles, is highly desirable 
particularly if we consider that in the not too distant future X- 
ray missions such as XEUS or Con-X will have, as one of their 
primary goals, the detailed characterization of relatively bright 
X-ray spectra. In some instances the advent of high resolution 
spectrometry will provide a valid solution. In others, where fea- 
tures that need to be characterized are broad, either because 
they are of non-thermal nature, such as the putative hard tails 
in clusters, or because they are smeared by other processes (i.e. 
the broad iron line observed in nearby A GN), solutions will 
have to be found elsewhere. Recent work by lDrake et al.l (l2007h 
provides an interesting starting point. 



Energy [keV] 



Fig. 4. Ratio of the ACIS-S3 (black) and pn (red) spectra accu- 
mulated in the l'-2' annulus over the best fitting MTP-gs-2Fe 
pn model. A re-normalization constant of 4.8% has been ap- 
plied to the ACIS-S3 spectrum. 



Technically this has been done by plotting the pn and the ACIS- 
S3 residuals with respect to the pn MTP-gs-2Fe best-fit model 
in the form of a ratio of the data over the model. A 5% re- 
normalization has been applied to the Chandra spectrum to 
facilitate the comparison with the pn spectrum. Interestingly 
the ACIS-S3 spectrum is significantly harder than the EPIC-pn 
spectrum. As shown in Fig. |4] while at 2 keV the two spectra 
have by construction the same intensity, at 5 keV the Chandra 
spectrum is in excess of the EPIC-pn spectrum by about 10%. 
The discrepancy is rather significant, considering that the non- 
thermal component intensity is of the same order. A compari- 
son of the pn spectrum with the MOSl and M0S2 spectra simi- 
lar to the one reported in Fig.|4]shows that in the 1-7 keV range 
both MOSl and M0S2 spectra are within a few % of the pn 
spectrum. We refrain from plotting the MOS spectra in Fig. |4] 
simply to avoid overcrowding. 

Since the background of the ACIS-S3 and EPIC cameras 
at energies below 5 keV is extremely small (see Fig. [TJ, the 
most likely origin for the difference is a cross-calibration is- 
sue on the effective areas in the 1-5 keV range. Interestingly 
cros s-calibration activ ities conducted by a number of workers 
(see lOavid et al.ll2007l) . mostly on cluster spectra, also point to 
a discrepancy between Chandra and XMM-Newton measure- 
ments. For hot clusters, (kT > 5keV) Chandra spec tra typically 



return higher temperatures than EPIC spectra (see lDavid et al. 



20071 figure on page 2). Further investigations have revealed 
that the discrepancy likely results from problems in the call 



bration of the Chandra high energy effective areas. iDavid et al 



(I2007h show that measuring the temperature of A2029 from 
the ratio of the He-like over H-like Ka lines and comparing 
the resulting continuum model with the observed spectrum, the 
latter is flatter than the former. More specifically, if the model 
is normalized to the observed spectrum at energies larger than 
5 keV, the latter is found to be about 8% below the model at 2 
keV. This discrepancy is remarkably similar to the one we find 
when comparing the pn and ACIS-S3 spectrum of Perseus (see 
Fig.H. 



5.1. EPIC estimate 

In this paper we have characterized the surface brightness of 
the non-thermal component in Perseus by resorting to a heuris- 
tic approach which can be divided into 2 steps: 1) inclusion 
of 2% systematic eiTors on spectra; 2) modification of spectral 
models. As far as the first step is concerned, we find a signifi- 
cant reduction of the surface brightness measured from the pn 
spectra when the 2% systematic error is included, no significant 
changes are found for the MOS measurements ( note however 
that a reduction is observed starting from a 3% systematic error, 
see Sect. 13.1.21 ). As far as the second step is concerned, the sta- 
bility of our results with respect to the inclusion of coiTections 
for some imperfections in redistribution matrix, effective area, 
energy scale and choice of astrophysical model indicate that 
our measures are relatively robust with respect to the residual 
systematic effects we have been able to identify. Comparison 
of MOS with pn results shows that the former favour a some- 
what larger value for the surface brightness of the non-thermal 
component than the latter, roughly speaking 2 against 1. The 
difference between the two measures is most likely caused by 
minor cross-calibration issues, possibly in the high energy re- 
sponse of either MOS or pn. Since we do not have firm evi- 
dence as to which of the two experiments is better calibrated, 
rather than privileging measures from one over the other, we 
take the common envelope of MOS and pn measurements (i.e. 
0- 5xlO"'^erg cm^^s^'arcsec"^) as our best estimate for the 
interval constraining the surface brightness of the non-thermal 
component. Some of our readers might find this range rather 
broad, particularly in light of the high statistical quality of our 
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data, we reiterate that the major source of indetermination here 
as elsewhere are systematic and not statistical errors. 

5.2. EPIC vs. ACIS 

Results from the analysis of EPIC spectra, either pn or MOS, 
cannot be reconciled with those obtained with ACIS-S3. As dis- 
cussed in Sect. |4] the difference is to be ascribed to a cross- 
calibration issue between EPIC and ACIS. Similar problems 
have been identified by a group of calibrators who have com- 
pared Chandra_and_2CMM-Newton observations of hot clus- 
ters ( David etal.l 20071). Recent efforts bv Chandra calibra- 
tors (IDavid et al.ll2007l : Ijerius et alj|2008l) have shown that the 
above difference follows from problems in the calibration of 
the Chandra high energy effective areas. 



5.3. Magnetic field estimates 

Under the assumption that the non-thermal emission originates 
from inverse Compton (IC) scattering of seed microwave and 
infrared photons by relativistic electrons, responsible of the 
synchrotron emission in the mini radio-halo, the Perseus core 
magnetic field can be estimated from the energy density of the 
seed photon field, Urad, the synchrotron and inverse-Compton 
luminosities, Lr and Lx, i.e. Ub - Urad ■ Lr/Lx, whe re 



Ufi is the energy density of the magnetic field. ISanders et al 



(120051) . from their measures of the non-thermal component, 
find magnetic field intensities ranging from several jjG at the 
very center, to ~ 1.0 juG at 8 kpc (0.40 an d ~ 0.1 yuG at 
40 kpc (2') (see Fig. 9 of ISanders et al.ll2005h . Our measure- 
ments can be used to provide a revised estimate of the B 
field. As discussed above we estimate the surface brightness 
of the non-thermal component to be somewhere between 
and 5xlO"'^erg cm^^s^'arcsec"^. We therefore start by assum- 



ing an upper limit of 5x10 erg cm s arcsec ', from this 
we subtract a value of 2.4xl0"'*erg cm^^s"' arcsec"^, as in 



Sanders et al.l (2005), to account for projection effects. The re- 
sulting surface brightness is used to derive lower limits for 
the magnetic field using e ssentially the same m ethod and ra- 
dio measures described in ISanders et alj (l2005b . In the 0.5'- 1' 
and r-2' annuli our X-ray estimate convert into lower limits 
of roughly 0.4 fxG and 0.3 fiG respectively. We note that these 
numbers have been derived through a series of rather drastic 
approximations, (i.e. we have a rather poor determination of 
the X-ray upper limit, the radio emission is measured at a fre- 
quency which is substantially larger than that at which non- 
thermal electrons responsible for the inverse Compton emis- 
sion emit via synchrotron, the B field is assumed to be con- 
stant, etc.) what really matter is that they are lower limits and 
not detections. 

Faraday rotation measures (RM) in cool cores on scales of 
tens of kpc are in the order of several hundreds to thousands 
(ITaylor et al.ll2002h . Estimates of magnetic fields from rotation 
measures have undergone some revision in the last few years 
with more recent estimates typical ly in the order of a few jjG 
( IClarkell2004t lEnBUn & Vo"gj|2006l) . In the case of the Perseus 
mini radio-halo Faraday rotation measures are available only 



on very small scales (ITaylor et al J 120061) . i.e. few tens of pc. 
RM estimates are in the order of ~7000 rad nT^ leading to 
B field values of ~ 25 yuG under the assumption the screen is 
localized in the ICM. This, however, appears to be unlikely as 
variations of 10% in the RM are observed on scales of ~ 1 pc 
(ITaylor et alJl20 06^. while ICM magnetic fields are ex pected to 



be ordered on significantly larger scales (few kpc : Tay lor et al 



20021: IVogt & EnBfinll2005l: lEnBlin & Vogtll2006l) . Application 
of the classical minimum-energy argument to the Perseus mini 
radio-halo data, leads to estimates for the central (i.e. r - ) 
magnetic field strength of ~ 7 fiG (IPfrommer & EnBlinll2004l) . 

IC esti mates of the magneti c field based on Chandra mea- 
surements (ISanders et alj|2005h are about an order of magni- 
tude below the RM and minimum-energy estimates detailed 
above. A similar discrepancy has been found when comparing 
magnetic field estimates on clu ster wide scales. Ther e IC mea- 
sures are in the 0.2-1 juG range (IRephaeli et al.l2008l and refer- 
ences therein), w hile RM estimates are about an order of mag - 
nitude larger (e.g. lCarilU & Tavloj2002l: iRephaeU et"ani2008l) . 
While the intricate nature of the astrophysical scenario lends 
itself to various possible explanations fo r the discrepancy (e.g. 
CarilU & Tayloill2002t fGovoni & Feretti..200 4). a more trivial 
alternative, namely that the IC estimates might be incorrect, 
should also be considered. There i s at least one object, Coma, 



where IC m agnetic field measures dFusco-Femiano et 



2004 , 120071) have been challenged dRossetti & Molendill2004 , 
20071) . Moreover first measures from the Suzaku mission on a 



laD 
ndi 



1999, 



few objects, (i.e. A3 667, Coma) are tu rning out to be lower lim- 
its on IC measures (IFukazawal 120071 and references therein). 
XMM-Newton measures presented in this paper seem to play a 
similar role in the determination of the Perseus core IC mag- 
netic field estimates. 



5.4. Future prospects 

It is unlikely that either Chandra or XMM-Newton will provide 
significantly better estimates of the intensity of the non-thermal 
component in Perseus than those reported here. Clearly experi- 
ments with sensitivities extending into the hard X-ray band are 
the most appropriate for these kind of studies. Recently the BAT 
experiment on board the Swift satellite has been used to deter- 
mine that an extrapolation of the non-thermal flux measured 
with Chandra to the 50-100 keV band overshoots the flux mea- 



sured with BAT by a factor of abo ut 4 CAiello et al.ll2008l) . The 
Suzaku high energy experiment (ITakahashi et al.l l2007b . lim- 
ited as it is by modest sp atial resolution and b y a non opti- 
mal background treatment (IKokubun et al.ll2007b . may provide 
some useful indications but will most likely not allow substan- 
tial advancement. Amongst missions that are currently under 
development Simbol-X (.Ferrando et al.ll2005h is arguably one 
of the most promising. The combination of large throughput in 
the 1-60 keV range, low instrumental background and EPIC- 
like spatial resolution, will allow a sensitive measurement of 
non-thermal components extending beyond the thermal cutoff. 
We have performed simulations of Simbol-X spectra based on 
our best fits for the l'-2' annulus trying out different values of 
the normalization of the non-thermal component to determine 
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how far down we might go. We started by taking a value of the 
surface brightness of 1 x lO^^^erg cm^^s"' arcsec"^. We find that 
the non-thermal component dominates thermal emission above 
^ 20 keV and remains above background emission up to ^ 30 
keV. Thus, if the non-thermal emission has an intensity of this 
magnitude or larger it will be detected rather easily. Repeating 
the same exercise with a surface brightness that is ten times 
smaller, i.e. IxlO^^^erg cm^^s^'arcsec"^, we find that formal 
fitting of the simulated spectrum still allows a detection of the 
non-thermal component at the ~ 3 - Act level, provided the ob- 
servation is longer than 3 x lO^s and cuiTent estimates of the 
Simbol-X instrumental background are within 10-20% of the 
real value. We note that in this case the relative intensity of the 
non-thermal component in the 20-30 keV band is about 10%. 
As we have learned from the analysis of our EPIC data, at these 
intensity levels, systematics, which are at this time unknown for 
Simbol-X, will play an important, possibly dominant role. In the 
specific case of Simbol-X measurements, the critical elements 
which need to be kept under control at the few percent level 
are: the effective area calibration; the cross-calibration between 
the high and low energy detectors and the background. These 
requirements, albeit challenging, are not beyond reach, partic- 
ularly if we consider that the formation flight strategy adopted 
by Simbol-X will entail some important advantages over pre- 
vious missions. Extensive calibration of the telescope effective 
areas can be performed in flight by comparing observations of 
calibration sources performed with and without the telescope in 
the optical path. Moreover direct illumination of the Simbol-X 
focal plane, with the telescope removed, will allow an in-flight 
verification of the detector quantum efficiency. 

Assuming that Simbol-X can indeed reach a sensitivity of 
IxlO^'^erg cm^^s^'arcsec"^ this will allow to detect IC com- 
ponents with associated B fields of roughly 1 juG, thereby pro- 
viding important constraints on the magnetic field in Perseus. 
Given the current magnetic field estimates fr om Faraday ro- 
tation measures on a few cool cores (few fiG: lEnBlin & Vogt 
200q) and the classical minimum-energy argument estimate on 
Perseus (~ 7^G: |Pfrommer & EnBlinl2004l) . even the exquisite 
sensitivity within reach for Simbol-X might be insufficient to 
detect the IC emission. 



6. Summary 

We have carried out a detailed analysis of a long EPIC observa- 
tion of the Perseus core in an attempt to detect and characterize 
a non-thermal component, our main findings may be summa- 
rized as follows; 

- systematic uncertainties play an important role in the char- 
acterization of the non-thermal component; in the absence 
of a strategy descending from first principles we have de- 
veloped a heuristic approach to include them in our analy- 
sis; 

- at variance with our preliminary estimates, we find that 
the non-thermal component is not detected; the sur- 
face brightness is determined to be smaller than ~ 
5xI0~^^erg cm^^s^'arcsec"^; 



our XMM-Newt on estimates ar e at var ia nce with Chandra 
estima tes from ISanders et alJ (120051) : ISanders & FabianI 
OOOTh and from our own analysis; the most likely expla- 
nation for the discrepancy between Chandra and XMM- 
Newton is a problem in the Chandra effective area calibra- 
tion; 

our EPIC based upper-limit on the surface brightness con- 
verts into IC magnetic field lower limits of ~0.4 /iG for 
the 0.5'-!' annulus and ~0.3 fiG for the l'-2' annulus; 
these measures are not in disagre ement with RM esti- 
mates on a few cool cores (few juG: lEnBIin & Vogtll2006h 
and the mi nimum energy estimate on Perseus (10 fiG: 



Pfrommer & EnBlin 2004); 



- in the not too distant future Simbol-X may allow detection 
of non-thermal components with intensities more than 10 
times smaller than those that can be measured with EPIC; 
nonetheless even the exquisite sensitivity within reach for 
Simbol-X might be insufficient to detect the IC emission 
from Perseus. 
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